THE NUMBER OF RHOMBUS TILINGS OF A SYMMETRIC 
HEXAGON WHICH CONTAIN A FIXED RHOMBUS ON THE 

SYMMETRY AXIS, I 

M. FULMEK AND C. KRATTENTHALER 

Abstract. We compute the number of rhombus tihngs of a hexagon with sides N, M, N, 
N, M, iV, which contain a fixed rhombus on the symmetry axis that cuts through the 
sides of length M. 



1. Introduction 

Let a, b and c be positive integers, and consider a hexagon of sides a, b, c, a, b, c whose 
angles are 120° (see Figure |l|.a). The subject of our interest is tihngs of this hexagon by 
rhombi of unit edge- length and angles of 60° and 120° (see Figure 0.b). (From now on, 
by a rhombus we always mean such a rhombus with side lengths 1 and angles of 60° and 
120°.) By a well-known bijection [^], the total number of rhombus tilings of this hexagon 
is equal to the number of all plane partitions contained in an a x 6 x c box. The latter 
enumeration was solved long ago by MacMahon Sec. 429, q ^ 1; proof in Sec. 494]. 
Therefore: 



The number of all rhombus tilings of a hexagon with sides a, b, c, a, b, c equals 

i+ j + k 
i + j + k 



i=l j=l k=l 

(The form of the expression is due to Macdonald.) 

A natural question to be asked is what the distribution of the rhombi in a random 
tiling is. On an asymptotic level, this question was answered by Cohn, Larsen and Propp 
[^. On the exact (enumerative) level, Ciucu and Krattenthaler and, independently, 
Helfgott and Gessel computed the number of all rhombus tilings of a hexagon with 
sides N, M, N, N, M, N which contain the central rhombus. (They were motivated by a 
problem posed by Propp [l^. Problem 1].) 

In this paper we solve an even more general problem, namely the enumeration of all 
rhombus tilings of a hexagon with side lengths A^, M, N, N, M, N which contain an arbi- 
trary fixed rhombus on the symmetry axis which cuts through the sides of length M (see 
Figure ^ for illustration, the fixed rhombus is shaded). Our results are the following. 

Theorem 1. Let m be a nonnegative integer and N be a positive integer. The number 
of rhombus tilings of a hexagon with sides N, 2m, N, N, 2m, N , which contain the l-th 
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A hexagon with sides N, M, N, A hexagon with sides N, M, N, 

N, M, N and fixed rhombus iV, M, N and fixed rhombus 

/, where N = 3,M = 2,l = l. I, where N = 3, M = 3, I = 2. 

The thick horizontal line indicates The thick horizontal line indicates 

the symmetry axis. the symmetry axis. 



Figure 2. 
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rhombus on the symmetry axis which cuts through the sides of length 2m, equals 



N N 2m 

X 

where the shifted factorial {a)k is defined by {a)k '■= a(a+l) ■ ■ ■ (a+k — l), k > 1, (a)o := 1. 



= 1 7 = 1 k = l 



Theorem 2. Let m and N be positive integers. The number of rhombus tilings of a 
hexagon with sides + 1, 2m — 1, + 1, + 1, 2m — 1, A^ + 1, which contain the l-th 
rhombus on the symmetry axis which cuts through the sides of length 2m — 1, equals 



r?^"') ^ \eJ{m + e){m + N-e){l-N) 

N+l N+12m-l . , . , , . 

^ n n n m^2- 

i=i j=i k=i 

In general, the sum in (1.2) and (1.3) (note that it is indeed exactly the same sum) 
does not simplify. (The only case where the sum is known to simplify is for A^ = 2n — 1, 
m = I = n, see |0, Corollary 3] and ||^, Theorem 12, Lemma 13].) Given that this is 
the natural "next" question is to ask what the quantities in (1.2) and (1.3) are 

"roughly" , or, more precisely, what the proportion of rhombus tilings, which contain the 
fixed rhombus, in the total number of rhombus tilings is asymptotically as the hexagon 
becomes large. Indeed, from Theorems |^ and 0, we are able to derive an "arcsine law" 
for this kind of enumeration. 

Theorem 3. Let a be any nonnegative real number, let h be a real number with < 6 < 1. 
For m ~ aN and I ~ bN, the proportion of rhombus tilings of a hexagon with sides 
N, 2m, A^, A^, 2m, A^ or A^ + 1, 2m - 1, A^ + 1, A^ + 1, 2m - 1, A^ + 1, which contain the 
l-th rhombus on the symmetry axis which cuts through the sides of length 2m, respectively 
2m — 1, m the total number of rhombus tilings is asymptotically 




;i-4) 



as N tends to infinity. 



This result is in accordance with the Cohn, Larsen and Propp result ||^, Theorem 1] 
that was mentioned above. The latter result does in fact give an asymptotic expression 
for this kind of enumeration for an arbitrary (semiregular) hexagon and an arbitrary fixed 
rhombus. (The result is even much stronger: The resulting function can even be used 
as a "density function"; consult the paper [Q for further information.) Indeed, the value 
'^7v,2m,7v(0, ^/3{—N/2 + /)), as defined in 0, Theorem 1], for m = aN and / = bN equals 

1 / a(l + a)-b(l-b) 

arccot 



vr \2^a{l + a)b{l -b) 
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which, in view of the formula 

o . f 1-2x2 
2 arcsm x = arccot , 

is exactly equal to ( |1.4|) , as it should be. 

In the next section we describe proofs of Theorems |I], ^ and ^ For the proofs of Theo- 
rems |1] and 1^ we build on the approach of . The basic ingredients are an application of 
Ciucu's Matchings Factorization Theorem [jl|, Theorem 1.2], the standard correspondence 
between rhombus tilings and nonintersecting lattice paths, and evaluations of the deter- 
minants which result from the nonintersecting lattice paths. As opposed to 0, here we 
have to take care of a few subtleties, which do not arise in ^ when the fixed rhombus is 
the central rhombus. Besides, we offer a few simplifications compared to 0. The proof 
of a crucial auxiliary lemma is deferred to Section 3. Finally, in Section 4, we point to 
further directions in this research, and hst a few open problems and conjectures. 

2. The proofs 

Proof of Theorems |1| and |2]. The proofs of both Theorems are very similar. We 
will mainly concentrate on the proof of Theorem 
There are four basic steps. 

Step 1. Application of the Matchings Factorization Theorem. First, rhombus tilings of 
the hexagon with sides A^, 2m, A^, A^, 2m, A^ can be interpreted as perfect matchings of the 
dual graph of the triangulated hexagon, i.e., the (bipartite) graph G{V, E), where the set of 
vertices V consists of the triangles of the hexagon's triangulation, and where two vertices 
are connected by an edge if the corresponding triangles are adjacent. Enumerating only 
those rhombus tilings which contain a fixed rhombus, under this translation amounts to 
enumerating only those perfect matchings which contain the edge corresponding to this 
rhombus, or, equivalently, we may consider just perfect matchings of the graph which 
results from G{V, E) by removing this edge. Clearly, since the fixed rhombus was located 
on the symmetry axis, this graph is symmetric. Hence, we may apply Ciucu's Matchings 
Factorization Theorem Theorem 1.2]. In general, this theorem says that the number of 
perfect matchings of a symmetric graph G equals a certain power of 2 times the number 
of perfect matchings of a graph (which is, roughly speaking, the "upper half of G) 
times a weighted count of perfect matchings of a graph G~ (which is, roughly speaking, the 
"lower half of G), in which the edges on the symmetry axis count with weight 1/2 only. 
Applied to our case, and retranslated into rhombus tilings, the Matchings Factorization 
Theorem implies the following: 

The number of rhombus tilings of a hexagon with sides N, 2m, A^, A^, 2m, A^ which con- 
tain the l-th rhombus on the symmetry axis which cuts through the sides of length 2m, 
equals 

2^-^R{S'{N,m))R{C{N,m,l)), (2.1) 

where S'{N, m) denotes the "upper half of our hexagon with the fixed rhombus removed 
(see the left half of Figure]^), where R{S'{m, n)) denotes the number of rhombus tilings of 
S'{m,n), where C{N,m,l) denotes the "lower half (again, see the left half of Figure ^), 
and where R{C{N,m,l)) denotes the weighted count of rhombus tilings of C{N,m,l) in 
which each of the top-most (horizontal) rhombi counts with weight 1/2. (Both, S'{N,m) 
and C{N, m, I) are roughly pentagonal. The notations S'{N, m) and C{N, m, I) stand for 
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C(2,2,2) 



The hexagons from Figure ^ cut along the symmetry axis, according 

to the Matchings Factorization Theorem. In the shaded regions, 

the "forced" (uniquely determined) rhombi are shown with thick lines. 

Figure 3. Hexagons, cut in two 

"simple part" and "complicated part", respectively, as it will turn out that the count 
R{S'{N,m)) will be rather straight-forward, while the count R{C{N,m,l)) will turn out 
be considerably harder.) 

It is immediately obvious, that along the left-most and right-most vertical strip of 
S'{N,m), the rhombi are uniquely determined. Hence, we may safely remove these strips 
(see the left half of Figure |^, the strips are shaded). Let us denote the resulting region by 
S{N — l,m). From ( p.l| ) we obtain that the number of rhombus tilings of a hexagon with 
sides N, 2m, N, N, 2m, N, which contain the l-th rhombus on the symmetry axis which 
cuts through the sides of length 2m, equals 

2^'^R{S{N - 1, m))R{C{N, m, I)). (2.2) 

Similarly, for the case of Theorem |^, we obtain that the number of rhombus tilings of 
a hexagon with sides N, 2m — 1,N,N, 2m — 1,N, which contain the l-th rhombus on the 
symmetry axis which cuts through the sides of length 2m — 1, equals 

2^-^R{S{N, m - l))R{C{N - 1, m, /)). (2.3) 

(See the right half of Figure |^. Note that in the case of Theorem ^, an applica- 
tion of the Matchings Factorization Theorem would directly give us 2'^~^R{S{N,m — 
l))R{C'{N,m,l)), with C'{N,m,l) the region as indicated in Figure ^. However, sim- 
ilarly to before, any rhombus tiling of the "complicated part" C'{N,m,l) is uniquely 
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Tilings for the "complicated parts" from Figure ^ interpreted as lattice paths. 
Figure 4. Lattice path interpretation 

determined in the left-most and right-most vertical strip of C'{N,m,l). Removing these 
strips then yields C{N — 1, m, I).) 

Step 2. From rhombus tilings to nonintersecting lattice paths. There is a standard 
translation from rhombus tilings to nonintersecting lattice paths. We apply it to our re- 
gions ^(A^, m) and C{N, m, I). Figure ^ illustrates this translation for the ( "complicated" ) 
lower parts in Figure |^. 

For the "simple" pentagonal part S{N, m) we obtain the following: The number 
R{S{N, m)) of rhombus tilings of S{N, m) equals the number of families (Pi, P2, • • • , Pn) 
of nonintersecting lattice paths consisting of horizontal unit steps in the positive direction 
and vertical unit steps in the negative direction, where Pi runs from {i,i) to {2i,i — m), 
t = l,2,...,N. 

Similarly, for the "complicated" pentagonal part C{N,m,l) we obtain: The weighted 
count R{C{N, m, I)) of rhombus tilings of C{N, m, I) equals the weighted count of families 
(Pi, P2, . . . , P/v) of nonintersecting lattice paths consisting of horizontal unit steps in the 
positive direction and vertical unit steps in the negative direction, where Pi runs from 
[2i — N — l,i + m) to [i, i) if i ^ I, while Pi runs from (2/ — N,l + m) to (/, I), with the 
additional twist that for i 7^ / path Pj has weight 1/2 if it starts with a horizontal step. 

Step 3. From nonintersecting lattice paths to determinants. Now, by using the main 
theorem on nonintersecting lattice paths [§, Corollary 2] (see also [l^. Theorem 1.2]), we 
may write R{S{N,m)) and R{C{N,m,l)) as determinants. Namely, we have 



R{S{N,m)) 



det 



m + 1 
m — i + j 



[2.4) 
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and 




Step 4- Determinant evaluations. Clearly, once we are able to evaluate the determi- 
nants in (2.4) and ( |2.5| ), Theorems |l| and ^ will immediately follow from (|2.2| ) and (|2.3|) , 
respectively, upon routine simplification. Indeed, for the determinant in (2.4) we have the 
following. 

Lemma 4. 

det (( \\ _-pr (Ar + m-i + l)!(i-l)!(2m + i + !),_! 

i<i!j<N\\m-i + j) J ~ f}^ (m + i-l)!(2A^-2i + l)! ' ^2 6) 

Proof. This determinant evaluation follows without difficulty from a determinant lemma 
in Lemma 2.2]. The corresponding computation is contained in the proof of Lemma 9 
in [§] (our determinant becomes the same as in when the order of rows and columns 
is reversed), and also in the proof of Theorem 5 in |^ (set r = N, Xg = m, B = 2, 
a + a — b = 2m there, and then reverse the order of rows and columns). □ 




On the other hand, the determinant in (|2.5| ) evaluates as follows. 
Lemma 5. 

- \y(m+i-j)\{N+j-2i)\ 

= ll (^ + ,_l),(2Ar-2i + l)! 11 + + 

(m)^^,nf=i(2j-l)! ^^ ^^JN\ (iV-2e)(i)e^ ^^.7) 



m YiW (202.-... ^ n e ; (rn + e) (.n + iV - e) (1 - iV) 

This determinant evaluation is much more complex than the determinant evaluation of 
Lemma ^, and, as such, is the most difficult part in our derivation of Theorems |l| and ^ 
We defer the proof of Lemma ^ to the next section. 

Altogether, Steps 1-4 establish Theorems ^ and |^. □ 

Proof of Theorem |3|. From MacMahon's formula ( |LTD for the total number of 
rhombus tilings together with Theorems |l| and ^ we deduce immediately that the propor- 
tion is indeed the same for both cases A^, 2m, A^, A^, 2m, A^ and A^ -|- 1, 2m — 1, A^ -|- 1, A^ -|- 
1, 2m — 1, A^ + 1, and that it is given by 

-mrr^)t^. (Ar-2e)(|) 



fnr-') hi ^ \e){m + e){m + N-e){\-N); 
We write the sum in (p.8|) in a hypergeometric fashion, to get 



(2Ar - 1)! {{m + l)N-if y^ {-N)e (1 - f )e {m), (-m - N), {\\ 

1)2/V-1 ^ 



(AT - 1)!2 (2m + l)2._i ^ (-f ), (1 - m - A^), (1 + m), (i - N), e! ' (2.9) 
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Next we shall apply Whipple's transformation (see |T8|, (2.4.1.1)]), which reads 

a, 1 + |, b, c, d, e, —n 



7-^6 



^,1 + a — b,l + a — c,l + a — d,l + a — e,l + a + 



(a + l)„(a — d — e + 1)^ 
{a - d + l)n{a - e + l)n 



4-^3 



a — b — c + l,d,e, —n 
6 + l,a — c+1, —a + d + e 



n 



;1 



Here, we used the standard hypergeometric notation 

(ai)fe 



s 



Cli, . . . , CL 

bi,...,h 



oo 

E 

k=0 



0'r)k k 



k\{bi)k---ibs) 



(2.10) 



(2.11) 



By setting a = —N, b = m,c = —m — N,d = 
as a limit case the following transformation: 

^ (-iV)e(l-f)e(m)e(-m-iV),(i) 



i e 

2' ^ 



-A^ + /,n = / - 1 in (2.10), we obtain 



e=ov f)e(l-m-iV)e(l + m)e(i-iV)ee! 

^ (-iV + l),_,(-/ + !),_, 
~ (-iV+i),_i(-/ + l),_i 



1 

-' 2 ' 



1 + m, 1 - m - A^, I' 



(2.12) 



Thus, expression 



turns into 



{2N - 1)! ((m + l)jv-i)' {-N + {-I + i)z_i 



■ 4-^3 



/ - A^, 1 - / 



2' 

1 + m, 1 — m 



AT 3' 



[N - 1)!2 (2m + l)2N-i i-N + i)^_i (-/ + l),_i 
Next we apply Bailey's transformation (see |1^, (4.3.5.1)]) between two balanced 4F3- 



series, 



4-^3 



a, b, c, —n 



e, /, l + a + 6 + c — e — / — n' 
4-^3 



(e - a)n{f - a)r 

(e)n(/)n 



X 



— ra, a,a + c — e — f — n + l,a + b — e — f — n + 1 
a + b + c — e — f — n + l,a — e — n + l,a~f — n + V 

with a = 1, 6 = 1/2, c = l — N,e = m + l,f = —m — N + l,n = I — 1. This gives 

(20! (2m)! {m + N- 1)! (m + A^)! (2A^ -21 + 2)! 



(2.13) 



4(/ + m - l)(m + - / + 1)(/ - 1)! /! (m - 1)! 

1 



X 



1 - /, 1, 1, 1 - / + AT 
2-l-m,2-l + m + N 



;1 



m! (AT - /)! (AT - / + 1)! (2m + 2N - 1)! [|, 2 - / - m, 2 - / + m + AT' ^^'^^^ 

for the ratio (|2.8| ). Now we substitute m ~ aA^ and / ~ bN and perform the limit A^ —>■ 00. 
With Stirling's formula we determine the limit for the quotient of factorials in front of 
the 4i^3-series in (2.14) as 2^a(a + 1)^6(1 -6)/(7r(a -b + l){a + b)). For the 4F3 -series 
itself, we may exchange limit and summation by uniform convergence: 



lim 4F3 

N^oo 



1-1,1,1,^-1 + N 
\,2-l-m,2-l + m + N 



;1 



1,1 

3 1 
2 



(1 - b)b 



A combination of these results and use of the identity (see 

arcsin \/z 



'a-b + l){a + b) 
I p. 463, (133)]) 



3 1 ^ 
2 



finish the proof. 



□ 



the number of rhombus tilings 
3. Proof of Lemma |5] 
The method that we use for this proof is also apphed successfully in [1^, |Tl|, |12|, |13 



(see in particular the tutorial description in [|T^, Sec. 2]). 

First of all, we take appropriate factors out of the determinant in (2.7). To be precise, 
we take 

{N + m~ i)\ 
{m + i-l)\{2N -2i + l)\ 
out of the i-th row of the determinant, i = 1,2, . . . , N. Thus we obtain 

AJ- {m + i- 1)! {2N -2i + 1)! 

X det /^/(^ + ^-i + l).-i(^ + ^'-2^ + 2k-.^^^^^ if^^A .3l^ 
i<M<^ ^(m + z-j + l),_i(iV + j-2^ + l);v-i+i iii = l) 

for the determinant in (2.7). Let us denote by D{m;N,l) the x A^-matrix underlying 
the determinant in (3.1), i.e., the (i,j)-entry of D{m;N,l) is given by 



D{m- N,l),j 



^m + i-j + l),_i(iV + J-21 + 2)^„ /+V'^' if ^ ^ ^ 



{m + i-j + l)j-i(A^ + j- 2i + l)N-j+i if i = I. ^^2) 

Comparison of (2.7) and (3.1) yields that (2.7) will be proved once we are able to 
establish the determinant evaluation 

LiV/2j 

det {D{m;N, I)) = JJ {{m + i)N_2i+iifn + i + l)N-2i) 



i=l 

^(.-.K.-.) jm)^^, n;=i(2j - 1)! ^ fN\ {N - 2e)(i)e 



NluWm^N-..,^ ^ ^ \eJim + e){m + N-e)C,-N), 

For the proof of (3.3) we proceed in several steps (see below). An outline is as follows. 
In the first step we show that nl=i ("^ + '^)N-2i+i is a factor of det {D{m;N,l)) as a 
polynomial in m. In the second step we show that nl=i^^("^ + i + \)N-2i is a factor of 
det {D{m; N,l)). In the third step we determine the maximal degree of det {D{m] N, /)) 
as a polynomial in m, which turns out to be (^^^) — 1. From a combination of these three 
steps we are forced to conclude that 

[N/2i 

det {D{m; N, I)) = JJ {{m + i)N^2i+i{m + t + l)N-2i) P{m- N,l), (3.4) 
1=1 

where P{m; N, I) is a polynomial in m of degree at most A^ — 1. Finally, in the fourth 
step, we evaluate P(m; A^, /) at m = 0, —1, . . . , — A^. Namely, for m = 0, —1, . . . , — [A^/2J 
we show that 

P{m-N,l) = {^-l)mN+{m^-m)/22{m^+m)/2-N+l^^^^ 

n7=~r(2j - 1)! nr=i(fe - +k-2m- i)\{^),_,{k - n)^_^ 
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Moreover, we show that P{m; N, I) = P{—N — m; N, I), which in combination with (3.5) 
gives the evaluation of P{m; N,l) at m = — \_N/2\ — 1, ... , —N + 1, — A^. Clearly, this 
determines a polynomial of maximal degree N — 1 uniquely. In fact, an explicit expression 
for P{m; N, I) can immediately be written down using Lagrange interpolation. As it turns 
out, the resulting expression for P{m; N, I) is exactly the second line of (3.3). In view of 
( ^.41 ), this would establish (3.3) and, hence, finish the proof of the Lemma. 

Before going into details of these steps, however, it is useful to record two auxiliary 
facts. The reader may, at this point, directly jump to Steps 1-4, and come back to the 
auxiliary facts when they are needed there. 

Auxiliary Fact I. There holds the symmetry 

det {D{m; N, I)) = det {D{m; N,N + l-l)). (3.6) 

This symmetry follows immediately from the combinatorial "origin" of the determinant. 
For, trivially, the number of rhombus tilings which contain the Z-th rhombus on the 
symmetry axis is the same as the number of rhombus tilings which contain the (A^ + 1 — /)- 
th rhombus. The manipulations that finally lead to the determinant det {D{m; N, I)) do 
not affect this symmetry, therefore det {D{m; N, I)) inherits the symmetry. 

This symmetry is very useful for our considerations, because for any claim that we want 
to prove (and which also obeys this symmetry) we may freely assume 1 < / < ["^^y^J or 
[^J <1<N, whatever is more convenient. 

Auxiliary Fact II. There holds the symmetry 

det {D{-N - m; N, I)) = (-l)^^^')-! det {D{m; N, I)) . (3.7) 

In order to estabhsh (13.71), we claim that D(m; N, I) ■ ((—1)-' (''■Z])) ■ ■ equals D(—N — 
m; N, I), except that all the entries in row / have opposite sign. Let us write R{N) for the 
matrix ((-l)-''(tJ))i<i^.<^- Since det{R{N)) = (-l)^^^^'), this would establish (g^). 

In order to establish this claim, we have to compute the {i, j)-entry in D{m; N, l)-R{N). 
For z 7^ Z, we have to show 



-!)■'■ X] (! !) + ^ - ^ + ^)k-i{N + k-2i + 2)j 
k=i ^ ^ 



= (_iV - m + 2 - J + l),_i(Ar + j-2t + 2)^._ /^^-"f-^'+^) , (3. 

and for i = I we have to show 

j 



{~iy I] _ J) + ^ - ^ + + k-2i + l)N-k+i 

k=l ^ ^ 

= _(_iV - m + ^ - J + l),_i(A^ -2t + j + 1)n-j+i- (3.9) 
Note that for j = 1, equation (3.8) collapses to 

_ (AT - 2z + 3)^_i^^^ = (AT - 2z + 3),v-i^^^^, 
which is of course true, so we may assume j > 1 in the following. 



THE NUMBER OF RHOMBUS TILINGS 
We convert the left-hand side of (3.8) into hypergeometric form, to obtain 



11 



iy-ii(_2m -N){N-2i + 3)^_i 3F2 



Next we apply the contiguous relation 



-2m - 1, 1 - j, 1 - i - 
-2m - N,3-2i + N 



m 



J (3.10) 



0,^1,^2 



3-^2 



a - 1, Ai, A2 



Bi, B2 



a,Ai + 1,^2 + 1 

5i + 1,^2 + 1 ' 



J (3.11) 

to the 3F2-series in ( p.lO|) . We want to apply the case where a = —2m — + 1 and 
Bi = —2m — A^. By inspection, in this case parameters cancel inside the two 3F2-series on 
the right-hand side of the contiguous relation, thus leaving two 2-Fi-series. So we obtain 



— i — m,l — j 
3-2i + N ' 



(j - 1)(1 - i-m) 
{-2m- N){N -2i + 3) 



2Fl 



2 — i — m,2 — j 
4:-2i + N 



;1 



(3.12) 

Each of the two 2-Fi-series can be evaluated by means of the Chu-Vandermonde-summation 
(see m, (1.7.7), Appendix (III.4)]), 

(c - a)„ 



a, —n_ 

1 -L 



{C)r 



(3.13) 



where n is a nonnegative integer. We have to apply the case where n = j — \ and 
n = j — 2, respectively. Since j > 1 in our case, the nonnegativity- condition is satisfied, 
and we obtain 

1 — 2m — A^, 1 — j, 1 — -j — m 



3-^2 



-2m - A^, 3 - 2i + A^ 



;1 



(-A^ + i - m - 1)(A^ + j + 2m - 1)(A^ - 2 + m + 2)j_2 



. (3.14) 



(-2m- A^)(A^-2i + 3)j_i 

Inserting this into ( p.lO| ) shows that (3.8) is true. 

In order to show (3.9), we convert the left-hand side into hypergeometric form, to obtain 

1 — z — m, 1 — j 



;_l)j(Ar-2^ + 2)^2i^i 



2 - 2i + A^ 



;1 



(3.15) 



Again, we can evaluate this 2-Fi-series by means of Chu-Vandermonde summation (|3.13|) . 
This proves (3.9). 

Now we are ready for heading into the details of Steps 1-4. 

Step 1. nl=i ("^ + Af-2i+i ^-5 a factor of det {D(m; N, I)). Here, for the first time, 
we make use of the symmetry ( p.6|) . It implies, that we may restrict ourselves to 1 < / < 

For i between 1 and [A^/2J let us consider row A^ — z + 1 of the matrix D{m; N, I). We 
see that the j-th entry in this row has the form 

^m + N-^-J + 2),_i {-N + 2z + j)^^^. ™p±i. 

Since (— A^+2z+j) = for j = 1, 2, . . . , N—2i, the first N—2i entries in this row vanish. 
Therefore (m + i)N-2i+i is a factor of each entry in row A^ — z + 1, i = 1, 2, . . . , LA^/2j. 
Hence, the complete product nl=i^^ ("^ + 'i)N-2i+i divides det {D{m; N, I)). 
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Step 2. nl=i^^("^ + i + ^)N-2i is a factor of D{m; N,l). Again we make use of the 
symmetry ( p.6|) , which allows us to restrict ourselves to 1 < / < ["^^y^J- 
We observe that the product can be rewritten as 

LiV/2j 

n + ^ + 1)^-2^ = nf=f + e + i)-Me,A^-e-l}. 
1=1 

Therefore, because of the other symmetry ( p.7|) , it suffices to prove that (m + e + 1/2)*^ 
divides D{m]N,l) for e = 1,2,..., \_N/2\ — 1. In order to do so, we claim that for 
each such e there are e linear combinations of the columns, which are themselves linearly 
independent, that vanish for m = —e — 1/2. More precisely, we claim that for k = 
1,2 e there holds 



J2( ■) ' (column (A^-2e + + of D(-e - 1/2; A^, /)) 

j=i vj/ 

(Note that this operation is really feasible for these values of e and k.) As is not very 
difficult to see (cf. [|1^, Sec. 2]) this would imply that (m + e + 1/2)^ divides D{m; N, I). 
Obviously, in order to prove (3.16) we have to show 

k 



^ [.j{e + i - j - k - ^ - N)M-2e+j+k X 



{2N + 3-2e-2i+j + k)2e-j-k-i ' ^ = (3.17) 
which is (3.16) restricted to the i-th row, 27^/, and 



fk\ 

Y^i .j(e + l + /-j-A;-| - N) M-2e-i+ j+ki'^N ~2e~2l+j + k + l)2e+i-j-k - 

(Ar_e-/-1), 

(_4)fc(iV_e-/ + l), 

X (e + 1 + / - J - i - N)N_2e-i+j{2N - 2e - 2/ + J + l)2e+i-j = (3.18) 

which is (3.16) restricted to the l-th row. 

Both equations (3.17) and (3.18) can be shown by the same kind of "hypergeometrics" 
(contiguous relation and Chu-Vandermonde) as was used for establishing ( ^.71 ). 

Step 3. det {D{m; N,l)) is a polynomial in m of maximal degree (^^^) — 1- Clearly, 
the degree in m of the (z, j)-entry in the determinant D{m; N, I) is j for i I, while it is 
j — 1 for i = I. Hence, in the defining expansion of the determinant, each term has degree 

(E,'i,j)-i = rn-i- 

Step 4- Evaluation of P{m] N, I) at m = 0, —1, . . . , — A^. Again, we make use of the 
symmetry ( p.6| ), and this time restrict ourselves to ["^^-^J < / < A^- On the other hand, 
by the symmetry ( p.7|) and by the definition (|3^) of P{m] N,l), we have P{m;N,l) = 
P{—N — m; A^, /). Therefore, it suffices to compute the evaluation of P(m; A^, Z) at m = 
0,-l,...,-LAr/2j. 
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What we would like to do is, for any e with < e < \N/2\, to set m = — e in 
compute det (-D(— e; A^, /)), and then express P{—e; N, I) as the ratio of det (-D(— e; N, I)) 
and the right-hand side product evaluated at m = — e. Unfortunately, this is typically a 
ratio 0/0 and, hence, undetermined. So, we have to first divide both sides of ( |3.4| ) by the 
appropriate power of (m + e), and only then set m = —e. 

Let e, < e < [N/2\ , be fixed. For A; = 0, 1, . . . , e - 1 we add 

^ M ■ (column {N +l-2e + k + i) of D{m; N, I)) (3.19) 

to column N + 1 — 2e + k of D{m; N, I). The effect is that then (m + e) is a factor of 
each entry in column N + 1 — 2e + k. This can be proven by exactly the same "hyperge- 
ometrics" as we used for the proofs of Auxiliary Fact II and Step 2. So, we take (m + e) 
out of each entry of column A^+1 — 2e + A;, /c = 0, l,...,e — 1 and denote the resulting 
matrix by Di{m; N, I, e). We obtain 



Di{m; N,l,e)ij := 



(m + z - J + l)j-i{N + j-2i + l)N-j+i if i = 



3.20) 



(3.21) 



if j < — 2e or j > — e, and 

Di{m;N,l,e)ij 

' {2e + i-k + m- N)N-2e+k{N -i + m + l)fcX 

(2A^ - 2e - 2z + 2A; + 3)2e-2fc-i if i ^ I 

{2N -2e-2i + 2k + 2)2e_2fc(2e + i-k + m- N)N-e~i+k>^ 

(m + e + l)j_e-i(l — i + m + N)^ ii i = I 

a N -2e + l < j < N - e. 

From what we did so far, it is straight-forward that we must have 

D{m; N, /) = (m + e)^Di(m; A^, /, e). 

A combination with ( |3.4|) gives that 

P{m;N,l) 

[N/2i 

= det{Di{m;N,l,e)) JJ ((m + i)e-i(m + e + l)^_i_e(m + i + i)^_2i)"^ (3.22) 

i=l 

Now, in this equation, we are able to set m = —e. Hence, in order to determine the 
evaluation of P{m;N,l) at m = — e, we need the evaluation of det (Di(m; A^, /, e)) at 
m = —e. 

Assuming I > ^^y^, we claim that the following is true: If e > A^ + 1 — /, then we have 



det {Di{-e;N,l,e)) = 0, 



(3.23) 
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N -2e N ■ 
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e + 1 


\o 
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- 

W + 1 - e 


N + 1- e 




1 


oV 

N +l-2e 



N - e 



N 



Figure 5. 



otherwise we have 
det (L>i(-e; A^,/,e)) 



N-e 



\{{k - 1)?{N + k-2e- - N)^.,. (3.24) 

i=i 

In order to estabhsh this we observe that Di{—e;N,l,e) has a block form which is 
sketched in Figure The figure has to be read according to the following convention: 
If a block is bounded from above by a horizontal line numbered ri at the left margin, is 
bounded from below by a horizontal line numbered r2 at the right margin, is bounded 
from the left by a vertical line numbered Ci at the bottom margin, and is bounded from 
the right by a vertical line numbered C2 at the top margin, then the block consists of the 
entries from rows i = ri, . . . , r2 and columns j = Ci, . . . , C2- This block form is easily 



established by routine verification directly from the definitions (|3.20|) and (3.21). 
Hence the determinant Di{m] N, I) factors as follows, 

det {Di{m; N, I)) = (-l)^(^-^) det(Q2) det(Qi) det(M). 



(3.25) 



Since Qi and Q2 are upper and lower triangular matrices, respectively, it is easy to 
evaluate det(Qi) and det(Q2)- We simply have to multiply all the entries on the main 
diagonal. We obtain 



UrJo k\{2e -2k- 1)\{N + k- 2e)\ 



det(Qi 



if / < AT + 1 
if / > + 1 



e, 
e. 



(3.26) 



Note that (|3]2|) follows immediately from ( lOSi ) and ( ^1261) . 

By assumption we have / > ^^^-^ > e, by assertion (p.23| ) (which is already proved) we 
may assume Z < — e, so we are sure to encounter row I in Q2. Multiplying the entries 



THE NUMBER OF RHOMBUS TILINGS 15 

on the main diagonal gives 

\ Af-2e+l N-2e 

det(Q2) = ( - I \{{j-l)\{N-2e-j + 1)n-,+i. (3.27) 

(Note that for e = 0, Q2 is the whole, unmodified matrix -D(0; A^, /).) 

For the evaluation of det(M) we employ Krattenthaler's Lemma (see IQ, Lemma 2.2]), 



det ((X, + AO ■ ■ ■ (X, + A,+i)(X, + B,) ■ ■ ■ (X, + B2, 



where Xi, . . . , X„, A2, . . . , An, and B2, . . . Bn are arbitrary indeterminates. 

The (i,j)-entry of matrix M is the {i,N — e + j)-entry of matrix D{—e;N,l). By 



assumption we have / > ^^^-^ > e, so we do not encounter row / in M. Hence, we must 



evaluate 

(^(zi^(Ar + 1-J + l)^_e+,_i(2X -e-2i+j + 2)e_,) . 



det 

By taking out {i — N)]\i_e from the i-th. row and (— 2)"^^-'^— |^ from the j-th column, we 
get 



Je^^((. + (-l-X + l))...(. + (- 



-X - 1 + ^)) 



X 



(2 + (-X-j + l))---(z + (-X-l)) , (3.30) 



which is precisely of the form required for (3.28) (set n = e, Xi = i, Aj = —1 — N + ^-^^ 
and Bj = —j — N + l). So we obtain 

ri(j - m^)j-i (3.31) 

for the (e x e)-determinant in (3.30). Multiplying this with the factors pulled out previ- 
ously, we have 

e 

det(M) = {-2Y-^-\e),\{{i - - N)^.,. (3.32) 

i=l 



By inserting ( p.26| ), ( |3.27| ) and ( p.32| ) into ( p.25| ) and simplifying the expression, we 



obtain (3.24). Thanks to (3.22), this establishes (3.5), and thus completes Step 4. 

This finishes the proof of Lemma ^. □ 

4. Open problems and conjectures 

In this paper, we computed the number of rhombus tilings of a hexagon with sides 
X, M, X, X, M, X, which contain a fixed rhombus on the symmetry axis. There are two 
questions which suggest themselves: 
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Question 1: What is the number of rhombus tihngs of a hexagon with arbitrary side 
lengths a, b, c, a, b, c which contain an arbitrary fixed rhombus? 

Question 2: What is the number of rhombus tihngs of a hexagon with side lengths 
a, b, c, a, b, c which contain several fixed rhombi? 

Whereas it is too much to expect "nice" answers to these questions in general, there 
is indeed hope for further nice results in special cases. We should mention that, on an 
asymptotic level, both questions are settled by the Cohn, Larsen and Propp result p. 
Theorem 1] that was mentioned in the Introduction. 

Regarding Question 1: Clearly, the approach that we used in this paper is not good 
enough for any generalizations in these directions. For, the use of the Matchings Fac- 
torization Theorem at the very beginning requires a reflective symmetry of the region 
that we are considering. So, under this approach we need to have a hexagon with sides 
A^, M, N, N, M, N with the fixed rhombus on the symmetry axis which cuts through the 
sides of length M. (The reader should observe that the Matchings Factorization Theorem 
cannot be used with respect to the "other" symmetry axis, i.e., the symmetry axis which 
runs parallel to the sides of length M. For, this symmetry axis does not have the required 
property that it "separates" the dual graph, cf. the statement of [|I], Theorem 1.2].) 

However, it seems that the approach taken by Helfgott and Gessel does allow to 
obtain further results in this direction. (Recall that they also obtained Theorems |I| and 
1^ for the cases where the fixed rhombus is the central rhombus.) In fact, a further result, 
using the approach by Helfgott and Gessel, has already been obtained by the authors and 
will be subject of a forthcoming article 0. There, we compute the number of rhombus 
tilings of a hexagon with sides N, M, N, N, M, N, the parameters N and M being of the 
same parity, which contain a rhombus which is "next to the center of the region". (Note, 
that, since and M have the same parity, there is no central rhombus.) 

Regarding Question 2, there are strong indications that, at least, there are "nice" results 
analogous to those of Theorems ^ and ^ if we consider rhombus tilings of a hexagon with 
sides A^, M, A^, A^, M, A^ which contain a given set of rhombi on the symmetry axis which 
cuts through the sides of length M. 

Suppose that we fix rhombi li,l2, ■■■ ,lr on the symmetry axis, and let L = {li,l2, ... ,lr}- 
If we want to know the number of all rhombus tilings which contain these rhombi, then 
we can use the same approach as we used for proving Theorems |] and |^. 

That is (see Section 2), we apply first the Matchings Factorization Theorem. It implies 
that the number of rhombus tilings of a hexagon with sides N, 2m, N, N, 2m, N which 
contain the rhombi from the set L equals 

2^-'-R{S{N - 1, m))R{C{N, m, L)), (4.1) 

where S{N — l,m) is the same "reduced" "upper region" as in Section 2, and where 
C{N,m,L) denotes the resulting "lower half, i.e., a region similar to C{N,m,l) (which 
appeared in Figure however, where along the "upper border" of the region the rhombi 
from the set L are removed. As in Section 2, the symbol R{C{N,m, L)) denotes the 
weighted count of rhombus tilings of C{N, m, L) in which each of the top-most (horizontal) 
rhombi counts with weight 1/2. A similar result (generalizing ( |2.3| )) holds for the case of 
a hexagon with sides A^, 2m — 1,N,N, 2m — 1,N. 
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Clearly, since R{S{N — 1, m)) is already known (see (2.4) and (2.6)), "all" we need is the 
weighted count R{C{N,m, L)). Again, the tiling problem can be translated into nonin- 
tersecting lattice paths, and from the lattice path interpretation we obtain a determinant 
for R{C{N,m, L)). To be precise, we have 

(N+m-iy. / N-j±l\ 




{m+i-jy.(N+j-2iy. LLl,>:^l.J (^^2) 

The reader should compare this determinant with the one in ( p.5|) . Now we do not 
have just one "exceptional" row, now we have r = \L\ "exceptional" rows. In fact, this 
determinant can be regarded as the "mixture" of two matrices, 

A=( {N + m-t)\ iv-,+1 - 

\{m + t-j)\{N + j -2t + iy} 2 . 

and 

{N + m- iy. 



B 



^m + i-3y{N + 3-2iyj,^^.^^- 

Now, to perform the task of evaluating the determinant, we would try to follow the 
proof of Lemma ^ in Section 3. Indeed, Auxiliary Facts I and II, and Steps 1-3 carry 
over, when suitably modified. In particular, the analogue of ( p.4|) reads: 

The determinant in ( [4.2|) equals 

lN/2\ lN/2\ 

JJ (m + i)N-2i+i n ("^ + ^ + \)N-2^ ■ Pim; N, L), (4.3) 

i=l i=r 

where P{m\ N, L) is a polynomial in m of degree at most r{N — r). 

The reader should compare with ( p.4| ) and observe the differences: First, the second 
product in ( [4.3|) is "smaller" than the corresponding product in ( |3.4]) . Second, and un- 
fortunately, the degree of P(m; N, L) is in general considerably larger than the degree of 
P(m; N, I). 

It is here where the problems start. Now it comes to carry over Step 4. At present, we 
do not know how to accomplish that. In order to determine P(m; A^, L) we would need 
r{N — r) + 1 evaluations of P{m; N, L). What we are able to do is to follow Step 4 in 
Section 3 and determine the value of P{m; N, L) at m = 0, —1, . . . , — A^. Unfortunately, 
for r > 1, this is not good enough. 

Still, we do believe that a reasonable formula for P(m; A^, L) can be found. In fact, it 
appears that, quite often, P(m; A^, L) does indeed factor further. At this point, we want 
to direct the reader's attention to the fact that this is also the case for P{m;N,l), the 
polynomial that was defined in (|3.4|) and which equals the second line of (3.3). For, from 
the expression given by the second hne of (3.3), it is immediately obvious that (m+/) n~21+i 
divides P{m; N,l). In particular, if / = 1 then P{m;N,l) factors completely into linear 
factors. 

A similar phenomenon seems to hold for P(m; A^, L) in general. Namely, for L = 
{1,2, ...,r} it appears that P{m;N,L) factors completely into linear factors (see Con- 
jecture |1] below). And, whenever we move a rhombus by 1 to the right then the degree of 
P(m; A^, L) in m increases by 2 (at least for large A^; compare Conjectures [^-^ below). 

We have worked out a few conjectures corresponding to "small" L, meaning that the 
numbers in L are small, which means that the fixed rhombi are far left in the hexagon. 



18 



M. FULMEK AND C. KRATTENTHALER 



Conjecture 1. Let m be a nonnegative integer and N and r be positive integers with 
N > r. The number of rhombus tilings of a hexagon with sides N, 2m, N, N, 2m, N, 
which contain rhombi l,2,...,r on the symmetry axis which cuts through the sides of 
length 2m, equals 

,i(.-i)(.-2iv) r^lf n ^'T\ {mi:iN-2i-iy\{m+i+i)^^,,^, 

f }L ^' M - 1)!! (m + i + l/2)^_,, 

2=1 J = l k = l •' 

The number of rhombus tilings of a hexagon with sides + 1, 2m — 1,A^+1,A^+1, 2m — 
1, N + 1, which contain rhombi 1,2, ... ,r on the symmetry axis which cuts through the 
sides of length 2m — 1, equals 

o|(.-i)(.-2iV) r^lf n ^?j imm-2z-l)\\{m + z + lU-2i-i 

M-r il (2z - 1)!! (m + z + l/2)^_,, 

N+lN+12m-l 

X n n n 'j^j^,^ 

i=i j=i k=i 

Conjecture 2. Let m be a nonnegative integer and N and r be positive integers with 
N > r + 1. The number of rhombus tilings of a hexagon with sides N, 2m, N, N, 2m, N, 
which contain rhombi l,2,...,r — l,r + l on the symmetry axis which cuts through the 
sides of length 2m, equals 

2i(r-i)(r-2JV) 3r(A^ - r) 

(2r - 1){2N -2r+l){m + r){m + N -r) 

/m+Af-l\2 N-2 
\ 2m / i=N—r 

^ (2z)!! {2N -2i~ 1)!! (m + ^ + l)iv-2.-i ^ + J + A: - 1 

1{ (2z-l)!!(m + z + l/2V_2. ^^^^ + ^' + ^-2' 

The number of rhombus tilings of a hexagon with sides N +1, 2m — 1, N + 1, N + 1, 2m — 
1, N + 1, which contain rhombi 1,2, . . . ,r — l,r + 1 on the symmetry axis which cuts 
through the sides of length 2m — 1, equals 



2i(r-l)(r-2iV) 



3r{N -r) 



(2r - 1){2N -2r + l){m + r){m + N-r) 

/m+Af-l\2 N-2 ^ 
V 2m ) i=N—r 

'j\ {2i)\\ {2N-2i-l)\\{m + i + ^r^ ^r^ + i + A; - 1 

11 (2i-l)!!(m + i + l/2V_2. fiM + + 
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Conjecture 3. Let m be a nonnegative integer and N and r be positive integers with 
N > r + 2. The number of rhombus tilings of a hexagon with sides N, 2m, N, N, 2m, N, 
which contain rhombi 1, 2, . . . , r — 1, r + 2 on the symmetry axis which cuts through the 
sides of length 2m, equals 

^\{r-l)ir-2N)^ (Qa (iV - r - 1)2 

64 (r - 1/2)2 {N -r- 1/2)2 {m + r)2{m + N - r - 1)2 
(m^ + 2Nm^ + [n^ + ^^^^ - 20r!±2i±5^ ^2 

^ ((20r-+i)jv-20r^-2r-5) ^^ + ^{2r - l)(2r + l)(2iV - 2r - 1)(27V - 2r + 1)^ 
rr^r j| 1 j^(2i)!!(27V-2i-l)!!(m + i + %_2.-i 



X 



rr-') ti (2> - 1)!! (m + > + l/2)„.,< 

i + i + k-2 



N N 2m J 

xnnnm^- («) 

i=i j=i k=i 

The number of rhombus tilings of a hexagon with sides + 1, 2m — 1,A^+1,A^+1, 2m — 
1, A'" + 1, which contain rhombi 1,2, . . . ,r — l,r + 2 on the symmetry axis which cuts 
through the sides of length 2m — 1, equals 

2i(.-i)(.-2iv)f5 (r)2 (AT - r - 1)2 

64 (r - 1/2)2 {N-r- 1/2)2 (m + r)2 (m + A^ - r - 1)2 



(rn^ + 2Nm^+{N^+^^^ 



20zf±2r+5\^2 



^ ((2o.+i)iv-2o.^-2.-5) ^^ + ±{2r - l)(2r + l)(2Ar - 2r - l)(2Ar - 2r + 1) 
rr^r j| 1 j=|(2z)!!(2Ar-2^-l)!!(m + . + l)^_2.-i 



2m+2N- 
2m 



) ^! (2^ - 1)!! (m + ^ + l/2)^_2. 

AT+l Af+1 2m-l 



X n n n jry^^- 

i=i j=i k=i ■' 

Without difficulty we could move on and work out further conjectures. Already these 
three conjectures do contain so many similarities, so that the "next" formula, the formula 
for fixing rhombi 1, 2, . . . , r — 1, r + 4, can, almost, be guessed right away. Still, these 
guesses do not seem to help for the determinant evaluations that would be needed to 
prove these enumerations. So we content ourselves with these three conjectures, and hope 
that they give enough evidence that in this area there are a lot of further beautiful results 
waiting to be unearthed, and proved, of course. 
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